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SUMMARY 


Self excited vibrations may cause severe limitation of the performance of 
turbomachines. Bearing forces, elastic hysteresis and forces from the fluid flow 
through clearances are known as origin. Forces from the leakage flow became im- 
portant with increasing performance because they increase with it. Theoretical 
approaches and extensive experiments were made to determine the dependence of the 
forces from the leakage losses and from rotating flow in the radial gaps. There- 
fore evaluation of the exciting forces for the investigated types of turbine 
stages should be possible. 


Former investigations of the vibration system set the mixed damping co- 
of the bearings to zero, the calculation model was partly extremly 
simplified. A comparison with a symmetrical rotor without negligence of the 
mixed damping coefficients shows the limits of the way to do so and presents 
all important parameters for stability. In certain respect the stiffness and 
damping coefficients itself are not important, only the relation of them to each 
other. With a given bearing and its coefficients as a function of the so-called 
Sommerf eld-number the relation between running and critical speed as well as the 
ratio bearing/rotor stiffness are the determining parameters for stability. Un- 
certainties in the bearing coefficients, which may result from differences bet- 
ween ideal test bearing and real rotor bearing are able to change the result in 
stability calculation in a wide range. 


INTRODUCTION 


6^cited rotor bending vibrations are of great importance for turbo- 
machinery. While it is possible to reduce vibrations caused by the unbalance of 
the rotor to a sufficient degree by balancing, the self-excited rotor whirl is 
only to be eliminated by design corrections. In thermal turbomachinery it 
seems there are three origins for self-excited rotor whirl. The first kind is 
the internal damping or elastic hysteresis, the second arises from the lubri- 
cating oil-film in the journal bearing (oil whip) and the third is the so called 
clearance excitation, induced by fluid flow in radial gaps. The two first causes 
depend on the running speed of the rotor, the vibrations appear quite suddenly 
at a certain speed limit. The clearance-excitation increases with increasing 
power of the turbomachine, therefore the vibrations appear suddenly at a certain 
threshold power. The frequency of these vibrations generally corresponds to the 
lowest critical speed of the rotor. The following considerations ignore the 
rarely important influence of elastic hysteresis. 
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VIBRATION EXCITING FORCES AT THE ROTOR 


Bearing- and clearance excitation forces are instationary external forces 
at the rotor. With xg and yg as the displacement in two rectangular coordinates 
the bearing force follows with 
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The spring- and damping coefficients c and d are defined as functions of the so 
called Sommerf eld-number. If the coupling values - that are the coefficients 
with mixed indices - are not zero, a displacement of the journal leads to a 
force rectangular to the displacement. That force is the cause for self-excited 
vibrations (oil whip) after reaching a certain running speed (ref. 1 and z; . 


The clearance excitation forces do not act on the journal but on the sur 
face of the rotor which is in contact with the working fluid. They result from 
the unsymme t r ical fluid flow through the radial clearances at rotor and blading 
which appears according to the eccentricity between rotor and casing. Thereby 
cross forces as well as moments are able to act on the rotor. Some theoretical 
and experimental investigations however show ” with respect to the actual pos 
sible accuracy of recording these effects - that it is sufficient to have only 
regard to the cross force as the essential factor in clearance excitation 
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where x and y are the coordinates of rotor displacement. 


( 2 ) 


Primary this effect was explained and calculated based on the leakage 
losses of blade tips or labyrinths. From eccentric po'sition in the casing non 
uniform tangential forces result at the turbine wheel, figure 1. The resultant 
of these forces - the clearance excitation force Q - acts rectangular to the 
displacement direction. In case of a compatible circumpolar vibration of the 
rotor the force Q anticipates 90° to the displacement. Is the loss of energy 
from damping less than the work done by the exciting force Q, then amplitudes 
will increase. At compressors Q acts in the opposite direction and is therefore 
able to excite countermoving vibrations. As tangential forces and power of the 
turbomachine are proportional vibrations begin at a certain threshold power 
(ref . 3 and 4) . 


Tests with shrouded turbine wheels showed considerable larger cross forces 
than explicable from leakage losses. The cause therefor is a non uniform pres- 
sure distribution above the shroudings as a result of the unsymmetr ical fluid 
flow through the sealing gap. Such cross forces arise at shaft glands too 
(ref. 5). By this means generally for the coefficient q in equation (2) must be 
put up 
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( 3 ) 


with qs as the part from clearance losses and q as the part from pressure dis- 
tribution. With U for the tangential force, p the pressure, the circumference 
angle and e the displacement of the shaft you get according to figure 1 

^ 2tt ^ b 2 tt 

Iq = ~ / COSU) dU; q = — J J sintf p r di/ dz (4) 

lf=0 ^ ® z=0 lf=0 ^ ^ 

where r is the radius of the shrouding, z the coordinate in axial direction and 
b the width of the shrouding, qg in equation (4) can be developed furthermore 
by known linear statements on clearance losses. Signifies = m • Ahfs/u the 

isentropic tangential force of the turbine stage with in the mass flow, Ah- the 
isentropic enthalpy difference and u the average circumferential speed you^can 
put up 


U. 

~ 1" (3) 

1 is the length of the rotor blade and K a factor which takes into account 
all essential parameters (ref. 4). For the part q several numerical procedures 
were investigated. You can put up (ref. 6) ^ 

^ ^Pb S (6) 

where s is the radial clearance, Ap the pressure difference at the shrouding 
and Kj^ a factor which depends essential — in case of given design of the seal- 
ing - only on the kinetic energy of the flow in circumferential direction 
divided by Ap^^ 

P c ^ 

r* - 

= 1a;— (7) 

D 

p is the density and c^^ the tangential component of the absolute outlet 
velocity of the stator blade flow. Figure 2 e.g. shows the theoretical func- 
tional relation between and C* for two shrouding designs. 

EXPERIMENTAL INVESTIGATION OF EXCITING FORCES 


To obtain data as reliable as possible, applicable for design of turbo- 
machines and to state theoretical assumptions, numerous tests at different tur- 
bine stage types (impulse and reaction type) were untertaken (ref. 7). Figure 3 
shows three investigated stage types. It was possible to measure the cross 
force Q directly dependent from the adjustable eccentricity, partly the pressure 
distribution above the shrouding too. Such "nearly-stat ic” force measurement 
was sometimes completed by real vibration tests. Measurements at shroudless 


305 



blading confirm, that the exciting force is calculable using the part q 3 . On 
the other hand tests with shroudings have shown the important influence of 
for that stage types. Because pressure distribution measurements were not pos- 
sible for all test series, the results were pointed out according to equation 
(5) as dimensionless "clearance-excitation-coefficient" 
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K„ of course combines the effects from leakage losses and pressure distribution. 
Figure 4 shows results of a 1-stage and 3-stage test rotor with 50 % reaction- 
blading and a stepped 3-chamber labyrinth at stator and rotor shroudings. K 2 is 
plotted over the load factor = 2Ah£g/u^. In figure 5 the results of a chamber- 
stage with low reaction is plotted. In both cases the part from leakage losses 
(corresponding Kg) is drawn in for comparison. 

The part of pressure distribution at the exciting force is particularly 
high for impulse- and low-reaction-stages because of the heavy swirl-flow in 
the gap between stator and rotor wheel. That shows figure 5 and is possible to 
deduce from figure 2. For that reason further tests were made measuring the 
pressure distribution too (ref. 8). The large influence of swirl flow suggests 
to reduce the exciting force by disturbing the circumferential component of 
flow. Therefore three labyrinth designs were investigated according to figure 6. 
For disturbing the swirl-flow axial sheet strips were mounted in chamber one 
(design B) and chamber one to three (design C) . In figure 7 the testing results 
are plotted. As shown it was possible to reduce the exciting force essentially 
by help of swirl-disturbing-sheets respectively to eliminate it. Anyway it is 
important to choose an optimal layout; obviously design B leads to considerable 
exciting forces with changed sign which may excite countermoving vibrations. 

For the rest at this stage the part of pressure distribution is exceeding; the 
part from leakage losses does not correspond with theory or is vectorial com- 
pensated - which is more probable - by the resultant force from pressure dis- 
tributions, in case B even over-compensated. 


THE STRUCTURE OF THE VIBRATION SYSTEM 


For investigation of rotor stability in the first instance a symmetrical 
laval-rotor with journal bearings may be suitable because of the large quantity 
of parameters. Thereby the following questions are to take into consideration: 

1. Simplification limits of the rotor model. 

2. Influences coming from the bearing properties, especially from data faults 
or manufacturing tolerances. 

Figure 8 presents the vibration system with mass m, rotor-stiffness Cj^; both 
bearings are equal. With equations (1) and (2) as well as the mass forces it is 
easy to put up the differential equation system for translational motion of the 
rotor-disk derived from the balance of the forces (ref. 8). 
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( 9 ) 


bo^' + b^y + h^y + b^x + b^y + b^x + b^y + b^x = 0 

The coefficients a and b are combinations of mass, exciting parameter as well 
as spring- and damping coefficients of the system. The usual starting solution 
leads to a characteristic equation of 6th degree for the complex natural fre- 
quencies which are imaginary at the stability threshold. By iterative solution 
the maximum of excitation is computed. It at least depends on the relation of 
stiffness rotor/bearing, to the type of bearing, the Sommerf eld-number and the 
relation between running speed and ciritical speed of the rigid supported rotor 
(tuning of system). 

A more simplified model (ref. 9) ignores the coupling coefficients of 
bearing damping . The components of bearing force coming from the coupling 
coefficients of stiffness are transposed to the centre of mass, the force from 
damping too. Stiffness of bearings and rotor shaft are considered as spring 
series. In that way we get the differential equations 

mx+dx+cx=-(c +a)v 
X X xy H/ 7 

•• j • /- ( 10 ) 

my + d^y + c^y = - (c^^ - q) x 

The spring coefficients c respectively c and damping coefficients d represent 
the effect of the parameters of rotor and bearings together. Coupling between 
X- and y-axis is pointed out by the coupling coefficients of bearings and 
clearance exciting force alone. For the stability threshold there is a closed 
solution. Also here the result only depends on relations between the parameters 
of the system. For sufficient stability it seems to be important - besides the 
damping - to choose the anisotropy of springs as great as possible. 

The more simplified system (10) compared with the exact system (9) shows 
the usability of the simplified system for qualitative conclusions; on the other 
hand - dependent on the given conditions - for quantitative investigations large 
differences may occur. Obvious the simplified calculation is only exact for an 
inelastic rotor. For real bearings the use of the solution according equation 
(9) is advisable because of the sensitivity of bearing coefficients as a func- 
tion of bearing geometry and manufacturing tolerance. Thereby all parameters 
can change their value the same time, and depending on bearing type and rotor 
stiffness this might result in very large changes of stability threshold. 
Unfortunately these changes are small for the not so interesting (because more 
unstable) cylindrical journal bearing while they are larger for other types of 
bearings. A clear prediction of the sign of the change is not possible. 

Exact calculation from equation (9) using coupling coefficients of bearing 
damping too gave following results: A larger stiffness ratio rotor /bearing 
leads to increase or decrease of stability threshold, increasing bearing load 
acts in the same unpredictable way. With increasing of the ratio running speed/ 
critical speed you will nearly always get a deterioration of the stability 
threshold. The omission of the coupling coefficients of damping influences the 
results in different amount and direction, therefore an omission of that bearing 
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influence is not reconunendable ; the fault can reach in same case a very large 
amount, see figure 9, wherein the ordinate S means the ratio of clearance 
excitation/rotor stiffness; S = q/c^^. 

The result s relat ing to the influence of Sommerf eld-number and stiffness 
ratio show that increasing bearing load can improve stability but it must not. 

In order to answer the question about the reliability of the calculation re- 
sults dependent on the bearing data systematically only one spring- or damping 
coefficient was changed by 10 % and the relative change of result was computed. 
Depending on the ratio of running speed/critical speedand ratio of stiffness for 
elliptical and three-wedge bearings we got differences up to 300 % to 400 % in 
stability threshold. These results can be amplified or moderated by taking into 
account all data faults, see table 1. A good prediction depends essentially on 
exact geometric bearing conditions - as used at the test bearings - that means 
exact manufacturing and mounting; on the other hand it is uncertain wether the 
exactness of bearing data is sufficient. Based on the knowledge of today the 
calculation looks pretty uncertain. But it is possible to use the tendencies 
of the results as criteria for design changes at the bearings, as was already 
stated by experience. 


CONCLUSIONS 


Clearance excitation forces originate from non uniform leakage losses and 
from non uniform pressure distribution above shroudings or shaft glands as 
function of the eccentricity of the rotor. Theroretical investigation of these 
effects was assisted by numerous experiments, so the clearance excitation 
coefficients should be reliable for the tested stage types. Swirl preventing 
sheets at shrouded bladings are able to reduce or eliminate the excitation 
forces . 

Considering the structure of the vibration system the great influence of 
journal bearings on the stability of the rotor is shown. For stability especi 
ally the anisotropy of stiffness is essential. The calculation model should be 
extended for multiple cased turbomachinery respectively shafts with more than 
two bearings, methods therefor are known. The strong dependence of the stabi- 
lity on bearing parameters leads to the idea to try to obtain the anisotropy 
not by help of the bearing geometry but by suitable design of the bearing 
pedestal . 
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Table 1. CHANGE OF MAX. SYSTEM EXCITATION (AS) IN PERCENT 
BY CHANGING ONE COEFFICIENT 10 % 

width/diameter W/D = 0.8 
ratio running/critical speed 
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FOR EXPLANATION OF THE CLEARANCE EXCITATION 



F igure 1 . 


CALCULATED COEFFICIENT OF CLEARANCE EXCITATION FORCE 
RESULTING FROM PRESSURE DISTRIBUTION 

a flat shrouding; b stepped shrouding with 3 sealing tips 



RELATED FLOW ENERGY C* 

E 


Figure 2. 
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CLEARANCE EXCITATION COEFFICIENT FOR STAGE B (50 % REACTION) 
axial clearance = 3.3 mm 

O, A stage displacement O 1 calculated 

□ rotor displacement • > from 

X stator displacement X leakage losses 


o 



LOAD FACTOR ¥ 


Figure 4. 
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CLEARANCE EXCITATION COEFFICIENT K 







Figure 6. 


TESTING RESULTS FOR 3 LABYRINTH DESIGNS (A,B,C) 

K 2 clearance excitation coefficient 
T load factor 

related flow energie 

□ from force measuring 

^ — integration of pressure distribution 

— — _ _ calculation from leakage losses 







SIMPLIFIED ROTOR MODEL 

(coupling springs and dampers not drawn) 



Figure 8. 


TEM EXCITATION FOR 4 BEARING TYPES 




